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Abstract 



Rare decay modes J /if) — > Dgir + + c.c, J/ip — > D~ir + + c.c. and J/ip — ► D°K° + 
c.c. are searched on a basis of 5.77xl0 7 J/ip collected with the BESII detector at 
BEPC. No signal above background is observed. We present upper limits on the 
branching fraction B(J/ip -► D s ir + ) < 1.4xl0 -4 , B(J/ip -► D--K+) <7.5xl0 -5 
and B(J/if; -► D°K°) < 1.7xl0 -4 with 90% confidence level. 



1 INTRODUCTION 



The hadronic decays and electromagnetic decays of J/ip have been widely 
studied, however weak decays of J/ip have not been studied in detail [1]. For 
the J/ip lying below the DD threshold, decays to DD are forbidden. Weak 
decays to one single charm meson accompanied by other non-charm mesons 
are allowed. The standard model predicts that the branching fractions of flavor 
changing processes via weak interactions, such as J/ip — > DX (where D stands 
for D s or D and X stands for ir or Kg) are at the level of 1CT 8 or below [2] 
and thus smaller than those of strong and electromagnetic decays. This is 
unobservable in current experiments. Extensions to the standard model, such 
as Top Color models [3], the minimal super-symmetric standard model with 
or without R-parity [4] and the two Higgs doublet model [5] , can enhance the 
branching fraction to ~1CT 5 . The study of J/ip — > D s ii + , J/ip — > D~n + and 
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J/ip — > D°K° provide an experimental check and serve as a probe to new 
physics [6,7]. Figure 1 shows the Feynman diagrams for these decay modes in 
the framework of the standard model, respectively. Charge conjugated states 
are implicitly included. In this paper, we'll present a search for J/ip — > Dgii + , 
J/ip — > D~7i + and J/ip — > DqK° in a sample of 5.77x10 events collected 
with the Beijing Spectrometer(BESII) [8] detector at the Beijing Electron- 
Positron Collider (BEPC) [9]. 




Fig. 1. Leading order electroweak Feynman diagrams for (a) J/tp — > D s tt + , (b) 
J/ip — ► D~ir + , and (c) J/i/j —> D°K° in the standard model. 



2 BESII Detector 



BES is a conventional solenoidal magnetic detector that is described in detail 
in Ref. [10]. BESII is the upgraded version of the BES detector [8]. A 12-layer 
Vertex Chamber (VC) surrounding the beryllium beam pipe provides track 
and trigger information. A forty-layer main drift chamber (MDC) located just 
outside the VC provides measurements of charged particle trajectories cover- 
ing 85% of 47r; it also provides ionization energy loss (dE/dx) measurements 
which are used for particle identification (PID). A momentum resolution of 
0.017a/1 + p' 2 (p in GeV/c) and a dE/dx resolution for hadronic tracks of 
~ 8% are obtained. Time of flight (TOF) of charged particles is measured 
with an array of 48 scintillation counters surrounding the MDC. The resolu- 
tion is about 200 ps for hadrons. Outside the TOF counters, a 12 radiation 
length, lead-gas barrel shower counter (BSC), measures energies and positions 
of electrons and photons. Solid angle is covered over 80% and resolutions of 
cte/E = 0.22/^/E (E in GeV), = 7.9 mrad, and a z = 2.3 cm are obtained. 
Outside the solenoidal coil, providing a 0.4 T magnetic field over the tracking 
volume, an instrumented flux return with three double-layer muon counters 
identify muons with momenta greater than 500 MeV/c. 
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Monte Carlo simulations are performed using a GEANT3 based program 
(SIMBES) with detailed consideration of the detector geometry and response. 
The consistency between data and Monte Carlo has been checked in many 
reactions in J ftp and tp(2S) decays with reasonable agreement. Details are 
described in Ref. [11]. 



3 Event Selection 

Due to a large J ftp hadronic decay background, non-leptonic decay modes 
do not offer ultimate sensitivity. In this analysis, D$ and D mesons are re- 
constructed via semileptonic decay modes: Dg — > <pe~v e , D~ — > K*°e~v e and 
D~ — > K°e~v e) D° — > K + e~v e . The neutrino is undetectable in the detector, 
but carries energy and momentum. Dg and D mesons can not be identified 
by their invariant mass. However, two body constraints can be applied in the 
mode J ftp — > DX. Thus they are identified using the momentum information 
of X meson. 

Four charged tracks are required in all selected modes, and the total charge 
must be equal to zero. In order to ensure well-measured momenta and reliable 
particle identification, all the tracks are required to be reconstructed in the 
main drift chamber with a good helix fit. Each track is required to satisfy the 
geometry selection criterion | cos 6 1 ! < 0.8 , where 9 is the polar angle, and 
must originate from the beam interaction region (except the decay daughters 
of Kg), which is defined by R xy < 2.0cm and \z\ < 20.0cm, where R xy and \z\ 
are the closest approach of the charged track in the xy plane and z direction. 

A kaon or pion candidate is required to satisfy Wx,n > 0.1%, where Wk,tt 
is the weighted likelihood of the hypothesis of kaon or pion, which combines 
the TOF and dE/dx information. To reduce the misidentification rate, the 
likelihood ratio Rr,-k = Wk,w/ (Wk + W n ) for kaon or pion is required to be 
greater than 0.7. The ratio of the energy deposit in the BSC to the momentum 
is used to construct the likelihood for electron identification, W e > 1% and 
Re = W e /(W n + W K + W-x) > 0.85 are required for electron candidate. The 
angle between the identified electron and the nearest charged track is required 
to be greater than 12° to remove backgrounds from photon conversion. Low 
momentum electrons and pions can not be unambiguously identified, leading 
to a momentum cut P e > 0.25 GeV/c for electrons. 

Neutral Kaons are identified via the decay to 7r + 7r~ . All pairs of oppositely 
charged tracks are assumed to be 7r + 7r~. The distance between the K s decay 
vertex and the beam axis is required to exceed 5 mm the in xy plane. Back- 
ground from J ftp's decaying to states with extra neutral particles is removed 
by constraining the number of isolated photons to be zero. An isolated photon 
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is a photon with an angle between the nearest charged track and the cluster 
of at least 22° and a difference between the angle of the cluster development 
direction in the BSC and the reconstructed photon emission direction of less 
than 60°. In addition the energy deposit in the BSC is required to be larger 
than 0.1 GeV. 

The neutrino of the semi-leptonic decay of the reaction J/ip — > DX remains 
undetected. A kinematic quantity U m i SS = E miss — P m iss is used to identify 
missing neutrinos, where E miss and P m i SS are the energy and momentum of 
the neutrino. Ideally U m i SS should be consistent with zero. A requirement 
of |£^rmss| < 0.1 GeV may remove backgrounds from J/ip decaying to K®, 
i] and partial n final states which can not be rejected through the initial 
selection. Events with missing energy due to misidentified pions are rejected 
by P m iss > 0.2 GeV/c. Missing particles are required to be in the sensitive 
region of the detector to further supress hadronic background. 

In the decay mode J/ip — > Dg7i + , Dg mesons are reconstructed through their 
decay to Dg — > 0e~z/ e , candidates are reconstructed from two oppositely 
charged kaons, namely: — > K + K~ . The invariant mass of candidates is 
required to be within 0.015 GeV/c 2 of the nominal mass. In the decay 
mode J/ip — > D~7i + , D~ mesons are reconstructed through their decays to 
D~ — > K* e~u e and D~ — > K°e~v e . The K*° candidates are formed out of K~ 
and 7r + candidates. Their invariant mass is required to fulfil \Mk-k — M K *o\ < 
0.060 GeV/c 2 . The K s candidates are formed out of n + and n~ candidates. 
The mass window is {M^^- — M K o\ < 0.020 GeV/c 2 . As far as the decay 
mode J/-0 — > D°K° is concerned, D° mesons are reconstructed through their 
decay to K + e~i , e . Figure 2 shows the invariant mass distribution of KK, 
Kir and nn from the various decay modes. 

Figure 3 shows the recoil momentum distribution of D$ or D meson for 
the three decay modes after all selection criteria. The two decay modes of 
D~ — > K *e~u e and D~ — > K°e~v e are combined in Fig. 3(b). 



4 Monte Carlo Simulation 



Detection efficiencies are determined by detailed Monte Carlo simulations. 
J ftp — > Dgir + , J/ip — > D~n + and J ftp — > D°K° production and decay, in- 
cluding the detector response are simulated. 50,000 Monte Carlo events are 
generated for each decay mode. Branching fractions of Dg, D~, and D° decays 
are taken from the world average values [12]. The detection efficiencies and 
branching fractions are listed in Table I. 
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Inv. Mass of K~7t + (GeV/c 2 ) Inv. Mass of n Jt + (GeV/c 2 ) 



Fig. 2. The invariant mass distribution of resonance candidates (a) <j) from 
J/ip -D5VT+, (b) from J/V> -> D°i?°, (c) if* from J/V> -» £>"7r + , 
D~ — > i^°*e~^ e , and (d)X^ from J/ip — > D~7r + , D~ — » lf e _ f e . Data is shown as 
dots with error bar, the expected signal shape from Monte Carlo simulated signal 
events is shown as a histogram. The mass cuts are illustrated by arrows and are 
discussed in the text. 

Table 1 



Detection efficiencies and branching fractions [12]. 



Decay mode 


Intermediate decay 


e 




J/iP -> D- s k+ 


D s cj)e-u e ,<j)^ K-K+ 


2.46% 


(2.4±0.4)% x (49.2±0.6)% 


J/ip -> D~7T + 


D~ -> K*°e-v e ,K*° -> 


1.90% 


(3.74±0.21)% 




D" -> K e~u e , K° s -> vr+vr- 


3.76% 


(8.6±0.5)% x 1/2 x (69.20±0.05)% 


J/^ -» -D°K° 


D° -> K + e-z^ e , Kg -> vr+vr- 


4.82% 


(3.51±0.11)% x 1/2 x (69.20±0.05)% 



5 Systematic Errors 



The systematic error of the branching fraction is dominated by uncertain- 
ties of MDC simulation (including systematic uncertainties of the tracking 
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(a) 





(b) 



(c) 



0.6 0.8 1 1.2 

Momentum (GeV/c) 

Fig. 3. Momentum distributions of recoil D$ or D meson, (a) J ftp — » D^tt + , 
(b) J/tp — > D~ir + , (c) J/*0 — > I? ^ . Data is shown as dots with error bars. 
The expected signal shape from Monte Carlo simulated signal events is shown as a 
histogram. 

efficiency and other requirements). The efficiency varies between 14.2% and 
19.9%. We estimate a systematic error of 2% for requiring the number of iso- 
lated photons to be equal to zero, using J/ip — > pn decays [13]. The system- 
atic uncertainty of false electron identification is estimated to be 5% for each 
electron. Differences in the pion and kaon identification between Data and 
Monte Carlo simulation indicate a systematic error of 1.5% for each track. 
The errors on the intermediate decay branching fractions of D$, D~ , D°, <f), 
and K*° are taken from the PDG [12]. The statistical error of the Monte 
Carlo sample is also taken into account. The total number of J/ip events is 
(57.7±2.7)xl0 6 [14], determined from inclusive 4-prong hadronic final states, 
and 4.7% is taken as a systematic uncertainty. All systematic errors are sum- 
marized in Table 2. 



6 Results 

Since no excess is observed for the reactions J/ip — > Dg7T + , J/ijj — > D~n + and 
J ' /ip — > D°K° is observed above background, upper limits for the three decays 
are calculated. Using a Bayesian method, the upper limits for the observed 
number of events at the 90% confidence level are 1.70 for J/ip — > Dgir + , 6.21 
for J/ip — > D~7i + and 4.61 for J/ip — > D°K° as shown in Fig.4. 
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Table 2 





J/V> — > DgTT + 


J/tjj -> D~7T + 


J/tp — > D"K U 




— > (f>e~v e 


D~ K* e-u<, 


D~ ->■ i^°e-^ e 


D° — > if+e~!/ e 


MDC Simulation 


15.6% 


14.2% 


19.9% 


15.5% 


PllOtOIl VCtrO 


2.0% 


2.0% 


2.0% 


2.0% 


e PID 


5.0% 


5.0% 


5.0% 


5.0% 


7T, if pid 


4.5% 


4.5% 


4.5% 


4.5% 


B(D S ,D) 


16.7% 


5.6% 


5.8% 


3.1% 


MC Statistics 


4.0% 


3.2% 


3.9% 


3.4% 


number of J/tp 


4.7% 


4.7% 


4.7% 


4.7% 


total 


24.7% 


17.7% 


22.7% 


18.3% 



The upper limits on branching fractions are calculated by 



B 



' up < N JH> eB(\ - 



a sys) 



(1) 



where riffl is the upper limit of the observed number of events at the 90% 
confidence level. Nj/^ is the total number of J/tp events, e and B are the 
detection efficiency and branching fraction, respectively. For the decay mode 
J/tp — > D~n + , sB is the sum of the products obtained from decay modes 
D~ — > K°*e~i> e and D~ — > K°e~i> e . a sys is the systematic error of the three 
decay modes, in the range of 18.3% to 24.7%. In the decay mode J/xjj — > 
D~7i + , D~ mesons are reconstructed from two decay modes, the systematic 
uncertainty is 21.2% weighted by their branching fractions. 

Table 3 

Numbers used in the calculation of upper limits on the branching fractions of J/ip — > 
D s tt + , J/V> -> D'tt + and J /j) -> D°K°. 





J/Y> - Dgir+ 


J/%l> -► D-TT+ 


J /iff -► L>°if° 


obs 
n UL 


1.70 


6.21 


4.61 


eB 


2.90xl0" 4 


1.83xl0 -3 


5.85xl(T 4 


Sys. Err. 


24.7% 


21.2% 


18.3% 


b(90%c.l.) 


< 1.4xl0 -4 


< 7.5xl(T 5 


< 1.7xl0 -4 



In summary, we have searched for the decays of J/ip — > D s n + , J /ip — > L> tt + 
and J/^ -> -D -^ using 5.77 x 10 7 J/^ events taken by the BESII detector 
at the BEPC e + e~ collider. No evidence for any of these decays is found. The 
final results for the 90% confidence level upper limit of the branching fractions 
are listed in Table 3. 

Acknowledgements 
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Event Number 



Fig. 4. Likelihood distributions for the observed number of events for (a) 
J/ip — > DgTi + , (b) J/ip — > D~it + , and (c) J/ip — > D°K°. The observed number of 
events at a Bayesian 90% confidence level for the three channels are indicated by 
arrows in the plots. 
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